Towards the Koch Snowflake Fractal Billiard: 
Computer Experiments and Mathematical Conjectures 



Michel L. Lapidus and Robert G. Niemeyer 

Abstract. In this paper, we attempt to define and understand the orbits of the Koch snowflake fractal 
billiard KS. This is a priori a very difficult problem because d(KS), the snowflake curve boundary of KS, is 
nowhere differentiablc, making it impossible to apply the usual law of reflection at any point of the boundary 
of the billiard table. Consequently, we view the prefractal billiards KS n (naturally approximating KS from 
the inside) as rational polygonal billiards and examine the corresponding flat surfaces of KS n , denoted 
by Sks„- I n order to develop a clearer picture of what may possibly be happening on the billiard KS, 
we simulate billiard trajectories on KS n (at first, for a fixed n > 0). Such computer experiments provide 
us with a wealth of questions and lead us to formulate conjectures about the existence and the geometric 
properties of periodic orbits of KS and detail a possible plan on how to prove such conjectures. 



1. Introduction 

The long-term goal of this work is to justify the existence and investigate, via computer-aided experiments 
and suitable mathematical arguments, a prototypical fractal billiard, namely, the Koch snowflake billiard. 

Since the Koch snowflake curve (the boundary of the Koch snowflake billiard table, see Fig. [T]) is nowhere 
differentiablc, it is certainly not clear from the outset that the Koch snowflake billiard — or, let alone, the 
associated billiard flow — is a legitimate mathematical object of study. Indeed, the first question that springs 
to mind is "How do you reflect a billiard ball off of a nondifferentiable boundary?" On the other hand, 
because the Koch curve is highly symmetric (indeed, it is self- similar) and can be thought of as some kind 
of 'infinite polygon' (see Fig. [TJ, it may not be entirely unreasonable to harbor such an expectation. 

In this paper, we provide experimental and geometric evidence towards this claim, as well as formulate 
conjectures and open problems concerning this challenging problem. In a nutshell, the Koch snowflake 
billiard is viewed as a suitable limit of (rational) polygonal billiard approximations; and analogously for the 
associated billiard (and geodesic) flows. We also identify a variety of (potential) periodic and quasiperiodic 
orbits of the Koch snowflake fractal billiard, and discuss some of their properties. 

In order to present our experimental results and conjectures about the snowflake billiard, we need to 
briefly recall several concepts from the theory of dynamical systems (more specifically, of rational polygonal 
billiards) and fractal geometry. 

A polygonal billiard B is given by a polygon P and the bounded region enclosed by P. For example, the 
square billiard is given by the square and the region bounded by the square. When studying a polygonal 
billiard, one is primarily concerned with the orbits of a pointmass moving at unit speed in B. We assume 
that our billiard ball experiences no friction and perfectly elastic collisions with the smooth portions of the 
boundary. Upon colliding with a smooth portion of the boundary, we reflect at the boundary according to 
Snell's Law, which says that the angle of incidence equals the angle of reflection; see, e.g., |GaStVo] . |Guj . 
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Figure 1 . Construction of the Koch snowflake from the equilateral triangle A (with sides 
of length 1); from left to right: KSo = A, KSi, KS2, KS3. The process continues ad 
infinitum, yielding the Koch snowflake curve, viewed here as the boundary of the Koch 
snowflake billiard table KS. For each integer n > 0, KS n denotes the nth prefractal 
polygonal approximation to KS. Note that KS n defines a rational billiard because its 
interior angles arc either 7r/3 or 47r/3. While this property is not required to be a rational 
billiard, the interior angles of KS n alternate between 7r/3 and 47r/3. (Depending on the 
context, KS„ also sometimes refers to the polygonal boundary of this rational billiard; the 
latter boundary is more correctly denoted by d(KS n ), however.) 

Chapter 9 of [KaHa , or [Ta . In general, if the billiard ball collides with a corner of the polygonal billiard 
table B, we must terminate the trajectory because we can not determine reflection at such points^] We call 
the map that describes the flow on the billiard P the billiard map. An important case of a polygonal billiard 
is a rational billiard and is defined to be a polygon P for which every interior angle is a rational multiple of 
tt; see, e.g., jMaTaj . [Sm] . [Vo] or [Zo] , 

Although a fractal is a shape that lacks many of the appealing qualities of polygons, so long as we are 
dealing with the right fractal, it is a shape that can be generated by polygons, a fact that is heavily exploited 
in this paper. The type of fractals we are interested in are self-similar fractal curves in R 2 that are nowhere 
diffcrcntiable and have infinite length. Specifically, the fractal billiard that we eventually want to define 
and investigate is the Koch snowflake billiard, denoted by KS; the boundary of the associated billiard table 
is the Koch snowflake curve, denoted by d(KS). The construction of the Koch snowflake fractal curve is 
given in Fig. [T] The prefractal KS n of KS is then an nth level polygonal approximation of KS. Like many 
fractals, the Koch snowflake curve d(KS) has a non-integer dimension. In the case of the snowflake, that 
dimension is log 3 4, which indicates that the curve d(KS) has infinite length. 

The Koch snowflake is not a rational polygon. The boundary of the Koch snowflake 'fractal billiard' is 
nondifferentiable, hence making it very difficult to properly define, and let alone analyze, the billiard map. 
However, the behavior of a pointmass may be anything but random, because, as was alluded to just above, 
KS has the very special property that there is a sequence of finite polygonal approximations {KS n }^L 
converging to KS such that for each finite n, KS n is a rational billiard; see Fig. [l] In order to overcome the 
limitations of the Koch snowflake, we may examine the prefractal billiard approximations KS n and attempt 
to make an argument in support of the existence of periodic orbits of the limiting fractal billiard KS, based 
on results concerning the rational polygonal billiards KS n . Therefore, in some sense, the snowflake curve is 
viewed as an (infinite) 'fractal rational polygon' and the associated billiard table KS as a 'fractal rational 
billiard'. 

Roughly speaking, a compact set F C R 2 is said to be self-similar if it can be written as a finite union of 
scaled (or rather, similar) copies of itself. More specifically, this means that F is a nonempty compact subset 
of R 2 and there exists a finite collection {Wj}i =1 of contractive similarity transformations of R 2 , with I > 2, 
such that F = U<=i Wi(F). It then follows that for any nonempty compact subset X of R 2 , we have that 
linifc^oo W k (X) — F (in the sense of the Hausdorff metric), where W(X) := ULi w *( x ) and W fc denotes 
the fcth iterate of W; see, e.g., Chapter 9 of |Faj . For this reason, the self-similar set F is called the attractor 
of the iterated function system {Wj}| =1 . 



We will see that there are polygons for which we can determine reflection at a corner, the equilateral triangle billiard 
being one. 
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Figure 2. Construction of the von Koch curve (left), often simply called the "Koch curve" 
here. At each stage of the construction, the middle third of each line segment is replaced 
by the other two sides of an equilateral triangle based on that segment. The self-similarity 
of the von Koch curve (right). The Koch curve is decomposed as the union of four pieces 
similar to the whole curve. 



Figure 3. The Koch snowflake curve is obtained by pasting together three copies of the 
von Koch curve (from Fig. [2]). 



The Koch curve is a perfect example of a self-similar set; see Fig. Strictly speaking, the Koch 



snowflake fractal is not self-similar; specifically, it is the result of pasting together three copies of the Koch 
curve; see Fig. [3] In light of this, we make a slight abuse of language by referring to KS (or rather, its 
boundary d(KS)) as a self-similar fractal. 

We can exploit the fact that the proposed billiard KS has a self-similar boundary. We notice that, in a 
sense, KS n is constructed from equilateral triangles. The billiard ball dynamics on the equilateral triangle 
billiard A = KSq are well understood. In particular, the periodic orbits of the equilateral triangle billiard 
A are classified in |BaUm| . The resulting theorem, combined with the fact that KS n can be embedded 
in a tiling Ta„ of the plane by equilateral triangles with side lengths 1/3™, will aid us in explaining our 
experimental results and formulating some of our conjectures. 

This paper is structured as follows. In Section 2, we discuss the necessary theory of rational billiards 
so that we may better present our results on the prefractal polygonal billiards KS n ; see Fig. [T] In §2.1| 
we discuss a useful tool for discerning periodic orbits of rational billiards, namely, the unfolding of billiard 
trajectories. In general, such an unfolding gives rise to a surface that depends on the orbit, but we also 
state known results for constructing a surface Sp which is independent of the periodic orbit of the rational 
billiard P. We recall in §2.2| that the flow on the rational billiard P corresponds to the geodesic flow on the 
corresponding surface Sp. In fact, as is explained in §2.3| vertices constitute singularities of the billiard map 
and certain copies of the vertices in Sp constitute singularities of the flow on the surface. We discuss what 
implications this has for the billiard ball dynamics. 



In \ 2.4 we view the prefractal polygonal billiard KS n as a rational billiard, and discuss the consequences 
of this fact for the billiard flow at the vertices of KS n . We also construct the corresponding surface for KSi, 
KS2 and KS3 and indicate the singularities of each surface. We show that, in agreement with the general 
theory of rational billiards, reflection at acute angles can be defined, but that reflection at obtuse angles 
cannot be determined unambiguously. 



In the last part of Section 2, namely [2.5 we discuss the results of [BaUm on the equilateral triangle 
billiard and their implications for the billiard KS n . In doing so, we explain how KS n can be embedded 
in a tiling Ta„ of the plane by equilateral triangles with side lengths 1/3™ (as was alluded to above) and 
give an equivalence relation on the collection of all periodic orbits of the rational billiard A. We do not 
seek to generalize the results of |BaUm| to KS n but instead use these results to provide us with a list of 
appropriately adjusted initial conditions for testing in KS n . 



2 See also [Fa, pp. xiv-xv] and, in a related context, [Lal,2], [LaNRGj . [LaPaj or [La-vF, §12.3.1 & §12.5]. 
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In Section 3, we present and discuss our experimental results on the prefractal billiard KS n by raising 
and sometimes answering a variety of research questions. In general, we examine the behavior of what we 
call induced orbits of KS n . An induced orbit is one for which the initial condition was appropriately adapted 
from an initial condition in the equilateral triangle billiard A. 

In |3.1[ we examine the behavior of periodic orbits of KS n in the collection 7„ of billiard paths determined 
by the initial conditions (x m id,ir/3), where x m id is a midpoint of a side of KS n and 7r/3 is the initial angle 
of the periodic orbit at the starting point x m id- We attempt to organize the orbits in j n according to 
their geometrical or dynamical behavior, and provide definitions for what we call Fagnano orbits, piecewise 
Fagnano orbits and primary piecewise Fagnano orbits of KS n . While these are very interesting orbits, the 
elements in j n do not constitute all of the induced orbits of KS n . 



In S3. 2 we examine the behavior of periodic orbits induced by the initial condition (xq,0q), where 
xq =/= x m id or 6q ^ 7r/3 and the orbit remains nonsingular, meaning that it does not hit any vertices of the 
billiard table. A particularly interesting example of such a periodic orbit is given by (x m id + Sx, 7r/6), where 
Sx is a suitable value (or vector) used to perturb the initial basepoint lying at the midpoint x m id of a side of 
KS n . Because of the nature of the equilateral triangle billiard, the reflection at the vertices of A = KSq can 
be determined. Consequently, in j |3.3[ these singular (and periodic) orbits are then used to induce singular 
orbits of KS n , for some n > 0. 

Finally, in i 3.4 we close Section 3 by discussing a seemingly uninteresting type of orbit, which we call a 
quasiperiodic orbit of KS n . Such an orbit can be considered as a 'rational approximation' to a given periodic 
orbit. For example, if (xq,9 ) is an initial condition of a periodic orbit, then a corresponding quasiperiodic 
orbit would have an initial condition (xi,a/b), where a/b is a rational approximation of 0q (obtained via a 
continued fraction expansion of 9q). Because the boundary of the billiard KS n is changing as n increases, 
the study of quasiperiodic orbits of KS n should enable us to define a suitable notion of 'quasiperiodic orbit' 
of KS. 

In Section 4, we conclude our paper by stating several conjectures based on our experimental results and 
provide a list of open problems and conjectures related to the proposed fractal billiard KS and a particular 
collection of its periodic orbits. In particular, we conjecture the existence of what we call the primary 
piecewise Fagnano orbit of KS, which we propose would be a suitable limit of primary piecewise Fagnano 
orbits of the prefractal approximations KS n . In the special case of pp^ , the 'primary piecewise Fagnano' 
orbit of Ksf\ we also conjecture that its 'footprint]^] on the boundary d(KS) is a self-similar subset of 
d(KS) C M 2 ; more specifically, it is the natural counterpart of the middle-third Cantor set strung around 
the Koch snowflake curve d(KS). 

The other conjectures and open problems stated in Section 4 arc clearly of a longer term nature. They 
concern, for example, the existence of the Koch snowflake billiard KS (as a proper mathematical object) and 
of the associated billiard flowj^jor (in the very long-term) the relationship between the length spectrum of 
the elusive Koch snowflake fractal billiard and the frequency spectrum of the corresponding Koch snowflake 
drum (e.g., [Lai 2 , LaNRG ,LaPa] , along with §12.3 and §12.5 of |La-vF| )[] 

In closing this introduction, we mention that the subject of billiards — and particularly, polygonal bil- 
liards or even, rational billiards — is an active area of research in the field of dynamical systems. Books, 
survey articles and research papers on various aspects of this topic include, respectively, [KaHa.Ta], 
[GaStVo,Gu,HuSc,MaTa,Sm,Vo,Zo] and [BaUm,GuJul 2,KaZe,Ma,Vel 3]. To our knowledge, none 
of these references deals with the elusive case of 'fractal billiards'. It is certainly true, however, that these 
works (and the many relevant references therein) have played a key role in motivating and guiding our 
investigations on the 'Koch snowflake billiard' (and its rational polygonal approximations) reported in the 
present paper. 



The presumed periodic orbit pp,^ of KS is induced by the Fagnano orbit of the original triangle A = KSo (i.e., the shortest 
periodic orbit of A), along with its appropriate counterpart in each billiard table KS n \ see j ]3.1| along with Conjectures |4.4| 
and 4.5 of Section 4. 

4 l.e., the subset of the boundary d{KS) consisting of all incidence points of the periodic orbit pp,!^ of KS. 
^as well as of the ass ocia ted ' fract al surface' Sks an d of t ne corresponding g eode sic flow (conjecturally equivalent to the 
billiard flow on KS); see f 2.2 and § 2.4 along with parts (i) and (iii) of Conjecture 4.7 
6 See Open Problem 4.8 



2. Flows on Rational Billiards and Corresponding Surfaces 



In this section, we want to discuss the properties of the billiard flow on KS n and describe how KS n 
can be viewed naturally as a rational billiard. A rational polygon is a polygon with interior angles that 
are rational multiples of tt. If the boundary of a polygonal billiard is a rational polygon, then we call the 
associated billiard a rational billiard. For us to discuss the flow on the associated closed surface, we need to 
define a few concepts Q 

Definition 2.1. If H is a regular N-gon, where N is a positive integer, then Dm is the group of symmetries 
of H and has cardinality 2N. It is called the dihedral group and is generated by the reflections in the lines 
through the origin that meet at angles ir/N. 

Definition 2.2. Let A(P) be the group of planar motions generated by the reflections in the sides of a 
polygon P. Furthermore, let G(P) denote the subgroup of the orthogonal group 0(2) consisting of the linear 
parts of the elements of A(P). 

When the polygon P has a connected boundary, then we can give an alternate characterization of it 
being rational. 

Proposition 2.3. Let P be a polygon with connected boundary. Then P is a rational billiard if and only if 
the group G(P) is finite. In that case, if the interior angles of the rational polygon P are written in the form 
irirrij/nj) for j = 1, ...,r, where rrij, nj are coprime positive integers and r is the number of vertices of P, 
then G(P) is isomorphic to the dihedral group Dn, where N is the least common multiple ofni,U2, ...,n r . 

If P is a rational billiard, then we have at our disposal a method for discerning periodic orbits. This 
method was first given in [KaZe] and reiterated in various forms in the literature; see, e.g., the survey 
articles [Ga StVo] , |Gu] , |MaTa] and |Zo] . Let xq be the initial position of the billiard ball and 8q be the 
initial direction of the billiard ball. Let x\ be the point on the side s\ of the billiard table at which the 
billiard ball collides, and 6\ be the angle at which the billiard ball reflects off of s\. If the orbit a is periodic, 
then we want to illustrate this by 'unfolding' the billiard table and its contents. 

2.1. Unfolding the billiard table B. To unfold the billiard table and its contents, reflect the contents 
of B in the side s±. This results in a segment continuing in the direction of 9q. Moreover, we see that we have 
created a situation where the trajectory passes through s\ and continues on to a segment that is collinear 
with the initial segment of the trajectory (see Fig. [4] for an example in the equilateral triangle billiard). 
Continuing this process, we see that we can unfold a billiard table and its contents. When the orbit is a 
collection of line segments, such an unfolding will be a straight line. Furthermore, the length of the unfolding 
is exactly the length of the orbit in B. When the orbit unfolds and terminates on an edge of a polygon P 
in the unfolding such that the terminal point is a copy of the initial point and the segment makes an angle 
with the side of P equal to the initial angle, then the original orbit of P is periodic. 

It is significant that we can unfold an orbit to determine the periodicity of the said orbit. In order to 
construct a closed surface out of the unfolding (which is determined by the periodic orbit) , we identify sides 
of the unfolding based on where the ball collides. In particular, we identify the terminal side with the initial 
side since the orientation of P and the original polygon P are the same. This results in a closed surface^] 
whose geodesic flow^ corresponds to the billiard flow on P. We can see this clearly in Fig. [4] However, this 
closed surface depends on the initial position and angle. We want to construct a surface S = Sp that is 
independent of these initial conditions. In addition, the geodesic flow on that surface Sp will turn out to be 
equivalent to the original billiard flow on P. 

2.2. The invariant surface Sp and its associated flow. Consider the product PxA(P). Essentially, 
P x A(P) contains 2N copies of our rational billiard table P. A graphical representation of P x A(P) can 
be given by the following. If we fixed a vertex of a polygon P such that reflection in the adjacent sides of 
P generates 2N copies of P, then the resulting shape is what is called a generalized polygon. If two sides of 

For an introduction to the theory of rational polygonal billiards from various points of view, we refer, e.g., to IGaStVol . 
[Gu| . [MaTkj . [Sm] . pTaj . fVo] or |Zo] . 

^By 'closed', we mean that the surface does not have any boundary, 
^which is actually a straight-line flow on the closed surface 
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Figure 4. Unfolding the equilateral triangle billiard: the orbit unfolds to a straight line 
segment (denoted by the dashed line). The terminal point on the upper right equilateral 
triangle in the unfolding corresponds to the initial point and the angles are the same. Hence, 
this corresponds to a periodic orbit of the equilateral triangle. 

the generalized polygon are translates of each other and/or the result of a rotation by 7r, then we consider 
these two sides to be equivalent. Then, upon modding out by this equivalence relation, we obtain a closed 
surface Sp := P x A(P) / ~ that does not depend on any particular orbit. 

Example 2.4. Consider the triangular billiard P with boundary given by the triangle with interior angles 
(37r/8, 7r/8, 7t/2). We reflect P in the sides emanating from the vertex with angle ir/8. Then, as is expected, 
we get a surface with 16 copies of the billiard P; indeed, N = lcm{8, 8, 2} = 8 and so 2N = 16. Therefore, 
the associated generalized polygon is an octagon. Two sides are equivalent if they are translates of each 
other. Modding out by this equivalence relation, we have obtained the associated flat surface Sp as shown 
in Fig. [5} Note, however, that since its genus is greater than 1, this surface cannot be properly embedded in 
the plane, but that for convenience, we represent Sp as though it could be done; see Fig. [5] 

The surface Sp is often referred to as a 'translation surface' or a 'flat surface' in the literature; see, e.g., 
|Ve3| . [HuScj, [MaTa], |Vo] and |Zoj . Indeed, the geodesic flow on Sp is nothing but the straight-line flow 
in the associated generalized polygon, with pairs of opposite sides identified as indicated above. The key 
fact concerning the geodesic flow on Sp is that it is (dynamically) equivalent to the billiard flow on P. 

2.3. Singularities of the billiard flow on P and of the geodesic flow on Sp. As mentioned 
before, the billiard map is only well defined at certain types of vertices of the rational polygonal billiard. We 
demonstrate under which conditions the billiard flow can be defined at vertices. If P is a polygon, then the 
corresponding surface contains copies of the vertices of P. These copies of the vertices are then considered 
singularities of the flow on the surface and are called conic singularities. 

In order to understand what a conic singularity is, we must first introduce the notion of a conic angle. 
A conic angle is an angle that measures the radians required to form a closed loop about the origin. In 
the plane, the conic angle is 2tt. We can form spaces in which the conic angle is not 2n. Suppose we had 
a space in which the negative upper half plane was glued to the negative lower half plane and the positive 
lower half plane was glued to the positive upper half plane of a second copy of R 2 . Suppose we make similar 
identifications but instead glue the positive lower half plane to the positive upper half plane of the original 
copy of M 2 . Then we have formed a space in which the conic angle about the origin is Air. Consequently, 
there are two types of conic singularities: removable and nonremovable. 

A conic singularity is called removable if its conic angle is 2-7T. The conic singularity is called nonremovable 
if its conic angle is 2m7r, for some integer m > 2. In the context of rational polygonal billiards, we determine 
the conic singularity of a vertex as follows. Geometrically, if v is a vertex of a rational billiard corresponding 
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Figure 5. The closed surface Sp corresponding to the triangle P(3ir/8,ir/8,ir/2); the as- 
sociated generalized polygon is an octagon. Pairs of appropriate sides should be identified 
as indicated, in order to obtain Sp, a surface of genus g — 2. Note, however, that strictly 
speaking, this planar representation is not valid. We illustrate such an embedding here as a 
heuristic device. 
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Figure 6. The vertex with angle measuring 2ir/3 in the triangle above has a conic angle of 
An. Consequently, the only way to recover a correctly oriented triangle by way of reflecting 
through sides adjacent to this vertex is to sweep out an angle of Air. 



to a nonremovable singularity, when we sweep out an angle of 2it by continually reflecting the rational billiard 
in a side, the orientation is not preserved. In Fig. [6] we see that a vertex with an associated conic angle of 
Air must sweep out an angle of Air before returning to a copy with the same orientation. We calculate the 
conic angle of a conic singularity as follows. If, for j — 1, 2, r, 7r(mj /rij) is an interior angle of the rational 



billiard P, as in Proposition 2.3 then the corresponding conic angle is 2mjir. 



Consequently, the geodesic flow on Sp does not see removable singularities and continues on unimpeded, 
but is immediately deterred from progressing across the surface when the geodesic intersects a nonremovable 
singularity. The geodesic flow on Sp cannot be determined at such a singularity, implying that the billiard 
flow on P at the associated vertex cannot be determined eitherFj 



This implication follows from the equivalence of the billiard flow on P and the geodesic flow on Sp. 
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Figure 7. The surface Sp corresponding to the triangle P(3ir/8, 7r/8, 7r/2); the correspond- 
ing generalized polygon is an octagon, shown on the right. Again, we commit the sin of 
suggesting that the surface can be embedded in the plane, this time to illustrate the am- 
biguity experienced by the billiard ball at nonremovable singularities. Once the opposite 
sides of the octagon have been pairwise identified, as shown in the figure on the right, one 
obtains the surface Sp. Note that periodic billiard orbits in the triangle (left) correspond 
to straight-line paths on the octagon (right), and hence to closed (i.e., periodic) geodesies 
on Sp. This figure demonstrates that the two orbits emanating from points near a vertex 
corresponding to a nonremovable singularity behave in a radically different way; the two 
orbits are shown in pink and orange, with the grey orbit representing what one would like 
to believe is the orbit emanating from the nonremovable singularity. 



Example 2.5. Consider again the triangle P(37r/8, 7r/8, 7r/2), as was done in Example |2.4| From the above 
formula, the conic angle of the conic singularity corresponding to the vertex with angle 37r/8 is 6n = 2 ■ 3tt. 
Consequently, this conic singularity is of the nonremovable type. This means that the flow on the surface Sp 
cannot be determined through the corresponding copy of the vertex on Sp. Moreover, we cannot logically 
define reflection at this vertex. The vertex at the center of our surface, however, is of the removable type; 
the corresponding conic angle is 2 • lw = 2ir; see Fig. [7] 

2.4. KS„ as a rational billiard. It was important that we developed an understanding of conic 
singularities, because, in some sense, the Koch snowflake is entirely comprised of singularities. Whether or 
not 'corners' of the Koch snowflak^] billiard can be shown to correspond to removable or nonremovable 
singularities of the flow remains to be determined. 

The prefractal billiard KS n is a rational billiard. Indeed, an interior angle of KS n is either equal to 
it/3 (acute) or 47r/3 (obtuse) (^] Consequently, N — lcm{3, 3} — 3 and so the corresponding surface contains 
2 • N — 6 copies of KS n . Furthermore, removable singularities of the surface correspond to acute angles 
of KS n and nonremovable singularities correspond to obtuse angles of KS n . This means that when a is a 
singular orbit of KS n , this orbit can be naturally extended to a periodic orbit if the singularities encountered 
in KS n were of the removable type; see Fig. [8] Likewise, if a is a singular orbit and the billiard ball encounters 
an obtuse angle of KS n , the flow on the billiard P cannot be naturally extended. 

In Fig. |J we have represented the flat surface SkSi ; an d m Fig. [IT] we have depicted the generalized 



polygons corresponding to SkSu Sks 2 an d Sks 3 - With sides properly identified, the latter generalized 



llr That is, vertices of the prefractal polygonal approximations KS n , for any n > 0. It is worth noting that the union of 
these vertices for all n > forms a countable dense subset of d(KS). 

12 Furthermore, note that except for n = (when KSq = A, the equilateral triangle), we always have both removable and 
nonremovable singularities (and in equal numbers). 




Figure 8. Acute angles have associated removable conic singularities. Reflection at acute 
vertices can then be defined. We illustrate here a singular orbit of KS\. 




Figure 9. The flat surface SkSi corresponding to KS\\ sides identified. This surface is 
independent of the flow and is the true surface, with sides properly identified. We cannot 
naively reflect in the adjacent sides of a fixed vertex and subsequently identify sides to 
produce the surface, because Sks„ nas genus g > 1 for all n > 0. There are six copies 
of the table KS\ in the associated flat surface, because 2 • lcm{3,3} = 6. Moreover, these 
six copies are generated by letting the dihedral group D3 act on KSi with an orientation 
designated by a vector that is not parallel to any of the sides; having such an orientation is 
a necessary condition for producing the correct number of copies of KS\. In Fig. |10| below, 
we illustrate the associated straight-line flow on Sks\ for a particular periodic orbit of KS\ 
induced by the Fagnano orbit J^o of KSq = A. 
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Figure 10. The unfolding of the periodic orbit of KS\ induced by the Fagnano orbit JF 
of KS Q = A. 




Figure 11. The generalized polygons corresponding to KSi, KS2 and If 63, respectively; 
sides not identified. These should be viewed as the first, second and third level approxima- 
tions of the limiting 'generalized fractal polygon' presumably associated with the billiard 
KS. When the proper identifications are made, such surfaces SkSi > Sks 2 an d $KS a can be 
viewed as the first, second and third level approximations of the limiting 'fractal flat surface' 



(as hypothesized in part (i) of Conjecture 4.7 1 



polygons become flat surfaces and can be viewed as the first three prefractal approximations of the 'fractal 
flat surface' (of infinite genus) Sksj which we will conjecture (in Section [4]) to exist as a proper mathematical 



object; see part (i) of Conjecture 4.7 



2.5. Equilateral triangle billiards and implications for KS n . The equilateral triangle billiard A 
is a very nice rational billiard. The billiard dynamics on A are well understood. Two orbits of an equilateral 
triangle billiard are said to be equivalent if each orbit has an unfolding that is a translate of the other and 
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the unique shortest closed billiard trajectory of A and is given by the equilateral triangle 
inscribed in A. The associated initial conditions are (x m id, 7r/3) , where x m id is the midpoint 
of one of the sides of A and 7r/3 is the corresponding initial angle of the path. 




Figure 13. Consider a tiling of the plane by the equilateral triangle A„ with side lengths 
1/3™. Denote this tiling by Ta„ . Then we can see that KS n C T\ n . This observation is 
illustrated here for the case when n = 1. 



each have the same length. Denote the equivalence class of an orbit a by [a]. In [BaUm , Baxter and Umble 
give a classification of the periodic orbits of A, up to this notion of equivalence. Roughly speaking, they 
partition periodic orbits into two collections, those periodic orbits with odd period and periodic orbits with 
even period^] Among orbits with even period, they then partition the orbits based on special criteria. If a 
is an even periodic orbit, then its equivalence class [a] has the cardinality of the continuum. The collection 
of odd period orbits is a collection of odd iterates of the Fagnano orbit J^o, {■^o k+1 } < k=o- R cca U that the 
Fagnano orbit JPq is defined to be the shortest periodic orbit of the equilateral triangle billiard A^]it is the 
unique such closed path and is represented in Fig. |12| Specifically, for every k > 0, there is only one element 
of the equivalence class [J?Q k+1 ], that element being ^Q fc+1 , the (2k + l)th repetition of 

Roughly speaking, the billiard KS n is comprised of many copies of A„, an equilateral triangle billiard 
with side lengths measuring 1/3™. More precisely, if Ta„ is a tiling of the plane by equilateral triangles with 



side lengths measuring 1/3™, then KS n can be embedded in Ta„; see Fig- 13 Because of this observation 
and the fact that all orbits of A can be unfolded in Ta„ (see Fi g. [l4| , it is reasonable to expect that periodic 
and singular orbits of A extend naturally to KS n . In fact, Fig. |14| illustrates that an unfolded periodic orbit 
of KSi is collinear with the unfolding of the Fagnano orbit J^o of A = KSq. 

3. Experimental Results 

We want to describe the 'anatomy' of the prefractal billiard KS n . For fixed n > 0, KS n is just a rational 
polygon. However, as n changes, the boundary of our prefractal billiard KS n changes. We want to be able 



^The period of a billiard orbit is the number of times which the pointmass collides with the boundary. Consequently, an 
even period orbit is one for which the number of collisions with the boundary is even, and similarly for an odd period orbit. 

^The name "Fagnano" for the shortest orbit comes from the name of the Italian mathematician Giovanni Fagnano whose 
work dates back to 1775. Indeed, Fagnano researched the existence of the shortest inscribed polygons of a fixed polygon P. 
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Figure 14. As noted in the caption of Fig. 13 KS n C Ta„- If we consider a periodic 
orbit "f n ,i in 7n as shown, then we see that j n j unfolds to a straight line. This unfolding is 
collinear with an unfolding of the Fagnano orbit J^o of A„. (See ^3.1 for the definition of 
7„.) As in Fig. 13 the situation is illustrated here for the case when n = 1. 



to talk more concretely about the parts of KS n that 'sprout up' and those that 'die off' as n increases. Parts 
of KS n that 'sprout up' are what we will call cells of KS n . The collection of segments of KS n -i that are 
removed to construct the cells of KS n will then be called the ghost of KS n -\. The ghost of KS n -i will not 
be present in KS n , but such language will help us discuss the behavior of a certain type of orbits. 

We want to stress that we will not be discussing orbits of KS n . Rather, we will be discussing what we 
call induced orbits of KS n . What we mean by 'induced' is that if (xo, $o) is an initial condition of an orbit of 
A, an induced initial condition will be (x' Q , 6 ), where Xq and x' Q are collinear in the direction of 9o and x' Q is 
a point on KS n . The new initial condition (xq,#o) is what we call an induced orbit of KS n . The reason for 
such language is that an orbit of K S n eventually induces an orbit of KS n +k for some k > 1. Such an orbit 
is induced when the billiard ball passes through a 'ghost' of KS n+ k-i into a cell of KS n+ k- We will discuss 
this in detail in §3.1[ but first we formally define what we have discussed thus far and give an illustration of 
these concepts in Fig. [15] 

Definition 3.1 (The anatomy of KS n ). 

(i) (Ghosts of KS n ). Let n > and {s n ,i}i=i be the collection of segments comprising the boundary 
d(KS n ) of the billiard KS n . Then, for 1 < i < 3 • 4™, the open middle third of the segment s n .i is 
denoted by g n ^ and is called the ghost of the segment s n ^. Moreover, the collection G n = {g n ,i\f=i 
is called the ghost set of KS n . The segments g n ^ are removed in order to generate KS n +\; see 



Figs.[T5];a)-(c). 

(ii) (Cell of KS n ). Let n > 1, 1 < i < 3 ■ 4"" 1 and G„_i be the ghost set of KS n -i. Consider the 
set G n -\ U KS n and the bounded region in G„_i U KS n given by an equilateral triangle with side 
lengths measuring 1/3™. Then this bounded region is what we call a cell of KS n . We denote a cell 



of KS n by C n .i\ see Fig. 15 ^d) (by definition, KSq — A has no cell). 




A 





(a) The ghost set of (b) The elements of the (c) Out of every side (d) Go U KS\ . The blue 

KSq = A, denoted by ghost set Go are re- there 'sprouts' two seg- arrows indicate the cells 

Go- moved. mcnts, giving rise to C\,i, 1 < i < 3 • 4° = 3, 

KSi. of KSx. 



Figure 15. An illustration of Definitions |3. 1 [ i) -(hi) in terms of KS = A and KS\. The 
ghost of the segment so,i, denoted by go^, is a middle-third segment of so,i and is removed 
from s ,i so that we may construct the cell of KS\. Then g oi is referred to as the ghost 
of the cell C Xi of KS l . 



(iii) (Ghost of a cell). Let n > 1, 1 < i < 3 • 4 n_1 and C n ,i be a cell of KS n . Then the ghost g n -i,i of 
the segment s n _i t j is called the g/iosi o/ i/ie ceZZ G ni j. In other words, the ghost of the cell C„,i is 
the middle-third segment g n -i,i of the segment s n -i,i of d(KS n -i) that is removed from s n _i,j in 



order to generate that portion of d(KS n ); see the caption of Fig. 15 



Example 3.2. If we consider the base of the equilateral triangle with sides having unit length, then the 



ghost of the base sq,i of A = KSq is given by go,i = (1/3, 2/3); see Fig. 15 a) 



Definition 3.3 (Compatible sequence of orbits). For each n > 0, let a n be a periodic orbit of KS n . Then, 
the sequence {cfc n }^Lo ' s sa ^°- to be compatible if each a n is induced by ao, the initial orbit of A. 

We performed computer simulations of induced orbits of the billiard KS n for n — 1, 2, 3. As a result, we 
now proceed to pose, discuss and sometimes answer various research questions regarding the induced orbits 
of KS n , for n > 0. 

3.1. 7„ Orbits of KS n . We define 7„ to be the collection of periodic orbits of KS n with an initial 
condition of the form (x m id, tt/3), where x m id is a midpoint of an arbitrary side of KS n . When 7r/3 is not a 
feasible direction, let 7r/3 be an angle relative to the initial side. Every orbit in j n can then be demonstrated 
to be an orbit induced by some periodic orbit of the equilateral triangle A; see Fig. [14] for the case when 
n = 1, and for an illustration of the orbits in ji , see Fig. |16| 

We want to know what constitutes a shortest orbit of KS n . Is the shortest orbit an element of 7„? If 
not, are all the orbits in "f n of the same length? Can we justify calling an orbit "Fagnano" in some sense 
that keeps with the 'shortest length' meaning of the name? We now proceed to answer these questions by 
examining our experimental results. Out of this will come a clear definition of Fagnano orbit of KS n and a 
better understanding of the behavior of the orbits in j n . 

Let 7i j € 7i. Consider the unfolding of 71,.; induced by the classic Fagnano orbit J^o of A" So = A, and 
consider the orbit a induced by the initial condition (x m id, f/6). Then, the orbit a is clearly shorter than 
the orbit 71^. Now, let us also consider the orbit 72J € 72 induced by ^o- Examine the unfolding of 71^ 
and the unfolding of an equivalent orbit that does not reside on the ghost of the initial side, but on some 
other one-third-segment of the initial side. Since equivalent orbits have the same lengths, we can see that 
this is also a periodic orbit of A that is shorter than 72 j- but still longer than a. From this, we deduce that 
no orbit in 72 will qualify as a shortest orbit of KS2- Consequently, for every n > 0, no element of j n is 
a shortest orbit of KS n . However, this does not preclude us from identifying the shortest orbits in 7„; see 

Fig. [HI 
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Figure 16. There are two orbits in 71. 




Figure 17. We see that as we construct KS2 from KS\, the orbit of KS\ induces an orbit 
of KS2- Then, constructing KS3 from KS2, the induced orbit induces an orbit of if S3. 
However, this need not always be the case. Certain orbits would remain unchanged as we 
construct KS n+ i from KS n for some n. Moreover, we can see that the trajectory in the 
first image passes through the midpoint of the ghost gk.i of a cell Ck+i,i of KSk+i for k < 1. 
This fact will be significant in Sj4j 




A 



/ \ 
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FIGURE 18. From left to right: A n KS , A n KS X , A n KS 2 , A n KS 3 . For each approx- 
imation KSk with k < 3, we are removing a subset of the ghost set Gk of KSk- The set 
A n KS is comprised of three ternary Cantor sets. 



Consider the subcollection of j n comprised of orbits with initial point x m id £ An KS n ; see Fig. 18 
In general, our results indicate that this subcollection is the collection of orbits with shortest length among 
all orbits in 7„. We denote this subcollection by & n and call it the collection of Fagnano orbits of K S n . 
Recall from our earlier discussion in §2.5| that jF consists of a single element, namely, the Fagnano orbit 



(also denoted &q); see Fig. 12 



We give a special name to the periodic orbit induced by J^o i n the prefractal billiard KS n (for each 
fixed n > 0). We call this particular orbit the primary piecewise Fagnano orbit of KS n (denoted by pp^ n ) 
because, for every fc < n, the initial segment enters a cell Ck,i of KSk (by passing through the ghost gu—i,i of 




Figure 19. In this figure, we see an example of a non-7„ periodic orbit of KSq, KSi, KS2 
and KS3, where we consider the same inducing condition (x m j<j) 0q) m each iteration of the 
prefractal billiards. For a fixed inducing condition (xo,#o)i other simulations (not shown) 
are suggesting that induced orbits of KS n are (qualitatively) quite convoluted. 



the cell Ck,i, as defined in Definition 3.1 'iii)) and, upon reflecting, subsequently forms a true Fagnano orbit 
in a cell Ck- In general, based on this terminology, we refer to an orbit in 7„ as a piecewise Fagnano orbit. 
In Section 4, we will often denote the collection of piecewise Fagnano orbits of K S n by p^ n (instead of j n ). 

3.2. Nonsingular non-7 n periodic orbits of KS n . Our characterization of periodic orbits of KS n 
in terms of what we have referred to as "7 n " and "non-7„" orbits is by no means a rigorous classification of 
the induced orbits of KS n . Our results indicate that elements in 7„ are not always induced by the Fagnano 
orbit of A and non-70 orbits of A do not always induce elements in j n ; see Fig. |19| Since non-70 orbits of 
A can induce j n orbits of KS n , is it possible that non-70 orbits of A also induce singular orbits of KS n l 
Can we demonstrate that there is some nonsingular non-7„ orbit that induces a nonsingular periodic orbit 
of KS n+k , for all k > 1? 

If n = and xq 7^ x m id, then oi( XOl7r /3), the periodic orbit of KS n with initial (or 'inducing') condition 
(xo,7r/3), is certainly a permissible element of the collection of non-7„ orbits of KS n . If, for some integers 
j,v > 1 such that < v < , we let xq = (z//3 J ,0) on the base of A, then Qtr XOt „/s\ induces a singular 
orbit of KS n for all n > j; see Fig. [20] Moreover, this singular orbit collides with an obtuse angle of KS n , 
meaning that we cannot make sense of the billiard flow with this particular initial condition. However, 
experimental results are indicating that for an initial condition (xq,9q) = (x m id,^/6), the resulting orbit (i) 



will be periodic and (ii) may be nonsingular, for all n > 0; see Fig. 21 

With regards to the billiard A, changing the basepoint of the initial condition (xo, #o) wm n °t affect the 
periodicity of the orbit, but may affect the length of the orbit. In our experiments regarding orbits in j n , 
we observed that changing the basepoint did not affect the periodicity of the orbit or the length of the orbit. 
This can be seen by unfolding the orbits. However, our experiments on non-7„ orbits indicated otherwise. 
One would like to see that orbits emanating from a side s n ^ with the same direction are equivalent, but 
this did not turn out to be the case. In particular, an example of a nonsingular non-7„ periodic orbit that 
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Figure 20. The billiard flow cannot be determined at obtuse angles. However, examina- 
tion of the corresponding geodesic flow on Sks„ for n = 1,2 allows us to draw 'possible' 
paths; note the ambiguity at obtuse angles. It is intriguing that the 'possible' choices that 
the billiard ball 'makes' at the triadic points {1/3, 2/3} and {1/9, 2/9} are exactly the same. 
This observation may allow us to force a logically consistent definition of reflection at non- 
removable singularities of the billiard map on KS n . In other words, this figure suggests 
that there may be a rigorous way of reducing the ambiguity the billiard ball experiences at 
obtuse angles of KS n . Because of the unique 'infinite' symmetry enjoyed by the snowflake 
boundary d{KS), we may be able to express the billiard flow on KS at nonremovable sin- 
gularities in terms of a (countably) infinite number of 'possible' paths the point mass may 
take. Again, we may be aided in this investigation by considering the (yet to be determined) 
geodesic flow on the associated 'fractal flat surface' Sks] see Conjecture |4.7| 

increased in length after a translation of the basepoint is one given by the inducing condition (x m id, t/6); 
see Fig. [21] Unlike the case of 7„-type orbits, if we shift the basepoint of the inducing condition by some 
small perturbation Sx, then the resulting orbit is longer in length; see Fig. [22] 

3.3. Singular orbits of KS n . Recall that a singular orbit of A is an orbit which collides with a vertex 
of the boundary dA. As we discussed in §2.3[ vertices of a billiard with associated conic angles measuring 
2ir radians are considered removable singularities of the geodesic flow on Sp. Since every vertex in A has an 
associated removable conic singularity, the billiard flow may be naturally defined at corners of A. 

Let a be an induced singular orbit of KS n . Since a is a singular orbit of KS n , a collides with cither an 
obtuse angle or an acute angle of KS n . Recall that acute angles of KS n constitute removable singularities 
of the billiard flow whereas obtuse angles of KS„ constitute nonremovable singularities of that flow. Conse- 
quently, we may not be able to always determine how a behaves at vertices of KS n . As n increases, KS n 
grows. Therefore, we may ask whether or not an induced singular orbit a remains singular as n increases. 
Because a can only be naturally extended at an acute angle vertex of KS n , is it possible to prove that 
induced singular orbits always collide with acute angles? Moreover, do singular orbits which only collide 
with acute angles at some stage n do so at every later stage n + k, k > 1? Even though our attention is 
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Figure 21. In this figure, we see an example of a non-7„ periodic orbit of KSo, KS\, KS2 
and AT S 3 . We consider the same inducing condition (x mid ,Tr/6) in each iteration of the 
prefractal billiards. 



focused on induced orbits, is it possible that there exist singular orbits of KS n that are not induced by any 
orbit of A? 

If a is an induced singular orbit of KS n , then experimental results indicate that a will always have to 
collide with a vertex of KS n +k, for all k > 0. An explanation for why this may be so is that KS n +k can 
be embedded in T&. k and the inducing initial condition gave rise to a singular orbit of A. However, our 
experimental results show that there is no consistency in which types of vertices are encountered by the 
singular orbit. In particular, it can happen that a singular orbit of KS n may only collide with acute angles, 
but the induced orbit of KS n +k will collide with an obtuse angle for some k > 1 and thus end there; see 



Fig. 23 for a concrete example. 

We note that because KS n can be embedded in Ta„ (see Fig. 13 1, it is reasonable to expect that every 
singular orbit of KS n is induced. 

3.4. Quasiperiodic orbits of KS n . When one considers a rational billiard P, the natural question to 
ask is " What are the periodic orbits of P?" One usually does not consider initial conditions that do not give 
rise to periodic orbits ™ Moreover, if one has identified periodic orbits of P, aperiodic "approximations" 
to these orbits may not be worth analyzing. Nevertheless, when the boundary is changing — as is the case 
in the recursive construction of the Koch snowflake billiard table KS via the rational polygonal billiard 
tables KS n — such approximations, however such a notion is defined, may be interesting. We call these 
approximate orbits quasiperiodic orbits. Recall that every orbit of a billiard is assumed to have a unit-speed 
parameterization. Formally, we define what we mean by a quasiperiodic orbit of KS n as follows. 

Definition 3.4 (Quasiperiodic orbits of KS n ). Fix T > 0, n £ N, and let denote the collection of all 
periodic orbits of KS n . An orbit q of the billiard KS n is called quasiperiodic if given any e > 0, there exists 



See, however, the well-known Veech dichotomy [Vel— 3] (as described, e.g., in [HuSc,MaTa,Vo]), according to which a 
direction is either periodic or uniquely ergodic in the rational billiard table. 
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Figure 22. In this figure, we see an example of what happens when we shift the initial point 
from x m id to x rn id + Sx, for some small perturbation Sx. We consider the same inducing 
condition {x m i,i + Sx,tt/6) in each iteration of the prefractal billiards. What is interesting 
in this particular example is that the billiard trajectory seems to want to enter into a new 
cell with each iteration of the prefractal billiard. This indicates that there may be a limiting 
object to consider, presumably, a billiard trajectory of KS. 

a finite partition of [0,T], to = < t\ < £2 < ■•• < £/ = T, and periodic orbits pj € & n such that 
\q(t) -Pj{t)\ < e, for all t 6 [tj-i,tj] and j € {1,2, ...,l}. 

As the definition would indicate, quasiperiodic orbits exist in any rational polygonal billiard table, but 
may have special significance in the Koch snowflake billiard. In Fig. |24[ we see an example of a quasiperiodic 
orbit of KS\ that is 'shadowing' the 71 orbits of KS\. In general, if we consider a rational approximation to 
the initial angle 8$, we obtain a quasiperiodic orbit of KS n . If n is changing and the billiard ball is passing 
through the ghost of a cell C n +k.i, then we want to have a way of discussing the 'stability' of the resulting 
induced quasiperiodic orbits of KS n +k, for all fc > 1. 
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Figure 23. The induced singular orbit only collides with acute angles of KS\. The same 
initial condition induces a singular orbit of KS2, whereas this time the orbit collides with 
obtuse angles of KS2- Insufficient numerical accuracy prevents us from investigating (via 
computer simulations) the induced singular orbit of K S3. 

Figure 24. The left image depicts the quasiperiodic orbit q. The next two images depict 
the behavior of q over two subintervals of time (t , ti) and (ti, t-i). In the middle image, the 
corresponding periodic orbit pi would be the one induced by the initial condition {x m u, n/3) 
and the corresponding periodic orbit pi would be induced by ((1/4, 0), 7r/3), with the base 
of the generating equilateral triangle A lying on the x-axis and the lower left vertex being 
located at the origin. 

Definition 3.5 (Stability of a sequence of quasiperiodic orbits). Let qo be a quasiperiodic orbit of A and po 
a periodic orbit shadowed by qo over some subinterval [tj—\,tj] C [0,T] (as in Definition 3.4 1 . Then we say 
that the sequence {q n }^Lo of quasiperiodic orbits of KS n induced by q is stable if the sequence {p n }^Lo of 
periodic orbits induced by po has the property that q n shadows p n for every n > 0p^| 

With regards to quasiperiodic orbits of KS n , there are two types of potential behaviors that we would 
like to investigate. Specifically, we are naturally led to ask the following two questions. Firstly, if we fix 
n > as well as a given periodic orbit p of KS n , is there a sequence of quasiperiodic orbits {qi}^ 1 of KS n 
that converges to p, and if so, in what sense? Secondly, for a fixed quasiperiodic orbit q of KS n , in what 
ways is the induced quasiperiodic orbit of KS n+ i qualitatively and quantitatively different from ql More 
precisely, if q n is a quasiperiodic orbit of KS n , n > 0, with a fixed rational angle a/b independent of n, 
and {<Zn+fc}fcLo ^ s a sequence of induced quasiperiodic orbits, is the latter sequence stable (in the sense of 
Definition 3.5)? 



Our experimental results thus far are indicating that we can answer all of these questions in the affir- 



mative. If {di/bi} < *L 1 is a sequence of rational approximations of 7r/3, then Fig. 25 describes what happens 



^Hcncc, in the sense of Definition 



3.3 



above, {q n }'^L is assumed to be a 'compatible' sequence of quasiperiodic orbits. 
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Figure 25. This figure depicts what happens as the approximate angle 8q becomes more 
and more accurate. The quasiperiodic orbit seems to dry up as the initial angle 8q — > w/3. 
Here, like elsewhere in the experimental results reported in |3.4[ the initial angle 6 is a 
rational approximation of the initial angle 9 of a given periodic orbit p, obtained via the 
continued fraction expansion of 8. 



as a,i/bi — > 7r/3. This suggests that for a fixed n, there is some notion of convergence of quasiperiodic orbits 
to a corresponding induced periodic orbit a of KS n . For example, in the present case of Fig. 25 (where 
n = 1), the corresponding sequence {qi}°Zi of quasiperiodic orbits KS\ seems to be converging (in some 
suitable sense) to the primary piecewise Fagnano orbit pp,^\ of KS\. As n increases, Fig. 26 illustrates that 
the quasiperiodic orbits break up over the boundary, but that the induced sequence of quasiperiodic orbits 
appears to remain stable, in the sense of Definition |3.5[ Furthermore, if we increase the time T, as given in 
Definition |3.4[ then Fig. [27] shows that the same behavior occurs. 



4. Conjectures and Open Problems 

We propose here several conjectures and open problems regarding the billiard KS. As we have done 
in the previous sections, we begin by first defining a few necessary terms. Our ability to investigate the 
proposed billiard KS has been limited by the finiteness of the computer. Hence, our conjectures deal mostly 
with "what happens in the limit". However, it is not exactly clear what we mean by limit. Hence, when 
we say 'limit', we are assuming a suitable notion of limit, most likely (depending on the context) an inverse 
limit of some inverse system or the Gromov-Hausdorff limit of a sequence of complete metric spaces; see, 
respectively, |HoYo| or [GrJ for an introduction to these notions. 

Once a suitable notion of limit has been defined, we want to consider the (admittedly, possibly empty) 
collection of periodic orbits of the billiard KS, which we denote by 3?. However, our focus thus far has been 
on induced orbits of KS n . We define ,f to be the collection of induced periodic orbits of KS, and phrase 
our conjectures in terms of elements in J* . As in Section [3] where we discussed the 'anatomy' of KS n , we 
now discuss the anatomy of the proposed billiard KS and the collection of induced orbits J 1 ' . We define the 
footprinf^ of an orbit to be the points of the orbit that belong to the boundary of the billiard KS. We call 
the ghost set of KS the collection G = IJ^Li the union of all ghost sets G n associated with KS n (see 
Definition pnfr)). 

Definition 4.1 (Footprint of an orbit). Let a be an orbit of a billiard P. Then, the footprint of a is the 
set of points of the orbit that belong to the boundary dP. 

Definition 4.2 (Ghost set of KS). If G n is the ghost set of KS n , then the ghost set of KS is defined by 
G = U^Lo Gn- 



This term is sometimes used in the literature on billiards; see, e.g., |BaUm) . 
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Figure 26. As the snowflake "grows", we see that the quasiperiodic orbit breaks up over 
the boundary, but does not devolve into an incoherent collection of segments. Rather, there 
is some sense of stability. That qualitative sense of stability is captured quantitatively in 
Definition 13.51 



Definition 4.3 (Self-similar orbit). Let a be a periodic orbit of the proposed billiard K S. Then, a is said 
to be a self-similar orbit if its footprint is a self-similar subset of d(KS) C R 2 (^] 

Note, in particular, that by construction, the sequence of primary piecewise Fagnano orbits {pp^ n }^=o 
is compatible, in the above sense (see Definition 3.3 1. 



Our first conjecture asserts the existence of a very special periodic orbit of the proposed billiard KS. 
We denote the (presumed) 'limit' of the orbits pp^ n by ppjF, and call this the primary piecewise Fagnano 
orbit of KS. We can see from Fig. 28 that there is an iterated function system (IFS)p*| denoted by F, 
which is producing scaled copies of pieces of the orbit, and such that for any n > 0, pp^ n — U"=o -^(^o)- 
Observe that for each finite n, the chronologjQ_of the orbit pp^ n can be easily recovered. Furthermore, if 
we consider the fixed point attractor of the IFSFj 



& = lim F n (J? ), 

n — >oo 

then this set & should correspond to the footprint of the proposed periodic orbit pp.'F of JC^PI Defining 
the orbit pp,^ to be is not very satisfactory, however, because such a definition does not provide useful 



^Here and thereafter, when talking about an IFS or a self-similar set, we are making the same slight abuse of language 
as in *|l]when referring to the Koch snowflake curve d(KS) as being 'self-similar', whereas strictlyspeaking, 8(KS) is obtained 
by fitting together three copies of a self-similar set, namely, the von Koch curve; see Figs. |2|and|3| 

^More specifically, in the spirit of footnote |l8[ this 'IFS' is given by three differe nt I FS's, denoted by Fx, F%, F3 and each 
comprised of two contractive similarity transformatio ns o f R 2 ; see the caption of Fig. |28| The associated attractor is therefore 
a 'self-similar set' of R 2 (also in the spirit of footnote 181. 



The chronology of an orbit is the order in which the billiard ball visits points of the boundary of the billiard table. 



See, e.g., Fa. for a detailed discussion on iterated function systems. 
22 Note that by construction, J? is a self-similar subset of d(KS) C M 2 . 
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Figure 27. Comparing this to Fig. 26 we see that for fixed n > 0, a quasiperiodic orbit of 
KS n fills more of KS n as T increases, where T is as stated in Definition |3.4| If we consider 
a sequence of induced quasiperiodic orbits {q n }^Lo, with q n a quasiperiodic of KS n , then 
this figure, along with Fig. 26 illustrates that the stability of {q n }^ = o (in the sense of 
Definition 3.5) is (qualitatively) independent of the time T. Note, however, that both the 
partitioning of [0, T] into subintcrvals, and the time intervals over which the shadowing of 
Pn by q n occurs, are in general dependent on n. 



information on how to recapture any sense of chronological order. Accordingly, we make the following formal 
conjecture. 

Conjecture 4.4. Under a suitable notion of limit, we conjecture that the chronology ofppJP can be naturally 
realized and hence, that ppJ^ can be considered as a true periodic orbit of the Koch snowfiake billiard KS. 

Significant analytical evidence in support of pp^f being a well-defined orbit of the billiard KS is the fact 
that (i) pp,j£ is a self-similar orbit (specifically, its footprint appears to be a self-similar Cantor set) and (ii) 
has finite length. 

Once we have established the existence of the periodic orbit pp,^ ', we expect that demonstrating the 
existence of other periodic orbits given as 'limits' of compatible sequences of induced periodic orbits {a n }^L 
may not be an unreasonable endeavor. Specifically, we conjecture that orbits emanating from the boundary 
of the proposed billiard KS at an angle of 7r/3 can be realized by a suitable generalization of the procedure 
(yet to be precisely determined) that would substantiate our claim that pp^ forms a periodic orbit of the 
Koch snowfiake billiard. 

Our next conjecture makes specific Conjecture |4.4| and the comments surrounding it. It also discusses 
some of the geometric properties of the presumed 'piecewise Fagnano orbits' of KS, including the 'primary 
piecewise Fagnano orbit' pp^ , the existence of which was asserted in Conjecture |4.4| 

Conjecture 4.5 (Existence of self-similar periodic orbits of KS). 

(i) (The primary piecewise Fagnano orbit of KS). The periodic orbit pp,^ of KS can be realized (and 
its chronology restored) as a suitable limit (possibly, an inverse limit) of the compatible sequence 
of primary piecewise Fagnano orbits pp^ n of KS n . Furthermore, the periodic orbit pp.'P is a 
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Figure 28. Motivation for building the primary Fagnano orbit of KS n by means of an 
iterated function system (IFS): the IFS produces the orbit pp,^ n of KS n by first contracting, 
rotating and translating the Fagnano orbit of A, so that what results are three copies of 
J^o appropriately appended to J^o- In actuality there are three IFS's (Fi, F2, F3) acting in 
unison to produce this picture. Specifically F\ acts on the upper left scaled copy of J^o to 
produce two scaled, rotated and translated copies of that piece of ppj^i. Likewise, F2 acts 
on the upper right copy of jFo to produce two scaled, rotated and translated copies of that 
piece of the orbit pp&\. Finally, F 3 behaves similarly on the bottom scaled copy of J^o- 
The union of all three images unioned with pp,^\ is then the orbit pp.^i ■ Continuing in this 
manner, we can produce pp^ n , for every n > 0. 



self-similar orbit (in the sense of Definition 4-.S); more specifically, its footprint is the natural 
self-similar middle-third Cantor set, strung around the boundary d(KS). 
(ii) (Piecewise Fagnano orbits of KS). More generally, an arbitrary 'piecewise Fagnano periodic orbit' 
a of KS can be defined and realized as follows. There exists a compatible sequence {a n }^=o °f 
periodic orbits such that for each n > 0, a n belongs to p^ n , and {a n }SJLo converges (in a suitable 
sense) to arj Furthermore, a is a self-similar orbit of KS. 



The following open problem is motivated by our earlier discussion in £3.2 and \ 3.3| 



Open Problem 4.6. On the base of A = KSq, does there exist x (other than the midpoint x m id) such that 
the sequence {a n }'^L of compatible non-j n periodic orbits induced by the initial condition (x Q ,ir/6) will (i) 
avoid all obtuse angle vertices of KS n , for every n>0 and (ii) converge to a periodic orbit a of KS? 

Moreover, we conjecture that if the initial angle #0 of the trajectory is not 7r/3, but such that (xq,9q) 
would have induced a nonsingular non-7„ periodic orbit a n of KS n for all n > 0, then the associated 
compatible sequence {a„}^L 'converges' to a periodic orbit a of KS; and correspondingly, the associated 
sequences of footprints of {a n }^ =0 converges to the footprint of a. In addition, we conjecture that the 



"^Recall from the end of f 3. 1 



that p^n stands for the collection of piecewise Fagnano periodic orbits of KS n 



^Furthermore, it has an obvious counterpart in the case where xq ^ x mi( ; is allowed to be on any of the sides of A. 
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footprint of a is a topological Cantor setp] A primary candidate for illustrating this conjecture is the 



periodic orbit induced by the initial condition (x m id + Sx, tt/6), for some small perturbation Sx; see Fig. 22 
A key step towards placing the above conjectures and open problems in a broader context and on firmer 
mathematical grounds must involve the theory of flat surfaces associated with rational billiards, along with 
its eventual extension to surfaces of infinite genus (called here 'fractal flat surfaces') presumably associated 
with fractal billiards)^] Ultimately, we are conjecturing the existence of a suitably defined billiard table KS. 
This amounts to establishing that (i) there is a corresponding flat surface Sks and (ii) the billiard flow on 
KS is equivalent to the billiard flow on Sks- I n particular, we conjecture that the 'limit' of the sequence of 
prefractal billiard tables {KSn}^^ would correspond to this billiard KS, and that similarly, the associated 
surface Sks would be the 'limit' of the sequence of prefractal flat surfaces {SksAuLo- ^ we can P rove 
that there are removable singularities in the proposed surface Sks, then it is not unreasonable to expect 
the conjecture to be true. Moreover, if the wonderful symmetry of the Koch snowflake curve d(KS) can 
be exploited so as to provide us with a way of dealing with the ambiguity the billiard ball experiences at 



non-removable singularities (see, e.g., the caption of Fig. 20 1, then we would be further justified in believing 
the conjecture to be a plausible statement. 

We summarize the main aspects of this central conjecture — and clearly, long-term open problem — in the 
following more specific form. Naturally, part of the difficulty in dealing with this problem will involve finding 
the appropriate notions of limits involved in the formulation of the conjecture^] 

Conjecture 4.7 (The fractal billiard KS and fractal surface Sks, along with the associated flows). 

(i) (The fractal flat surface Sks, along with the geodesic flow). The sequence {Sks„}'^' = o of prefractal 
flat surfaces associated with the rational billiard KS n converges (in the Gromov-Hausdorff sense, 
see |Grj ) to a surface of infinite genus Sks, called the 'fractal flat surface' associated with KS. 
Correspondingly, the 'geodesic flow' on Sks can be realized as a suitable limit (possibly, an inverse 
limit) of the geodesic flows on the surfaces SKS n - 

(ii) (The Koch snowflake fractal billiard KS, along with its billiard flow). The fractal billiard KS can 
be defined as a suitable limit (possibly, an inverse limit) of the prefractal polygonal rational billiards 
KS n : essentially^^ this means that the billiard flow on KS can be realized as a suitable limit (also 
possibly, an inverse limit) of the corresponding billiard flows on the billiard tables KS n . 

(iii) (Geodesic vs. billiard flow). Finally, the geodesic flow on the fractal flat surface Sks is equivalent 
to (and, at first, may be used as a suitable substitute for) the billiard flow on the Koch snowflake 
billiard KS. 

We close this paper by stating the following very long-term problem, which is directly motivated by 
the questions raised in |La2] concerning the relationship between 'fractal drums' and 'fractal billiards', 
including the Koch snowflake drum and billiard. See, in particular, [La2j , Conjecture 6, p. 198, itself 
motivated by Conjectures 2 and 3, pp. 159 and 163-164, respectively; see also §12.5.3 of [La-vF . For 
information regarding the Koch snowflake drum and other 'fractal drums' (viewed as 'drums with fractal 
boundary'), see, e.g., [Lai 2,LaNRG,LaPa], [La-vF,§12.3 & §12.5], and the relevant references therein. 
For information regarding trace formulas (including the Gutzwiller and Chazarain trace formulas) in various 
contexts, we refer, e.g., to [Gzl,2], [Chj . |Co| and |DuGnj . 

Open Problem 4.8 (Fractal Billiard vs. Fractal Drum.). Once the existence of the Koch snowflake billiard 



KS has been firmly established (as hypothesized in Conjecture can one formulate, and eventually estab- 
lish, a suitable 'fractal trace formula' in this context? Presumably, the latter would be a fractal counterpart 
of the Gutzwiller, Chazarain and the Selberg trace formulas in this context, connecting the length spectrum 



2 ^I.e., it is a perfect and totally disconnected subset of the boundary d(KS). 

2 ^The flat surface Sp of a rational polygonal billiard P always has a finite genus. Moreover, a simple calculation based on 
the known formula for g(Sp) (see, e.g., Lemma 1.2, p. 1022 of MaTa:) shows that g n = g(S xs„ ) — > oo as n — > oo. Hence, it 



1.7 



is natural to expect that the presumed limiting surface Sks nas infinite genus; see Conjecture 

2 ^For example, even if the notion of 'inverse limit' is suitable for formulating part (ii) or the end of part (i) of the conjecture, 
one would still need to specify the maps involved in the definition of the underlying inverse systems; see, e.g., |HoYo| . 

2 ^It is clearly true (and well known) that viewed as a sequence of compact subsets of M 2 , the billiard tables KS n converge 
in the sense of the Hausdorff metric (or, equivalently, in the Gromov-Hausdorff sense) to the Koch snowflake billiard table KS. 
Hence, the real issue concerns here the associated billiard flows. 
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of the snowflake billiarc^\ and the eigenvalue (or frequency) spectrum of the corresponding snowflake fractal 
drum. Moreover, can one address the same problem for other fractal billiards (once they have been properly 
defined) and the associated fractal drums? 

We hope that the study of such open problems and conjectures, aided by a suitable combination of 
computer experiments and theoretical investigations, will enable us in the future to better understand the 
elusive nature of the Koch snowflake billiard and, eventually, of a variety of other fractal billiards. 



Added note. For a version of the paper with crisper pictures, please email the second author. Such a 
file will, however, be considerably larger and require between thirty and sixty minutes to download over a 
standard 56Kbps connection. 
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